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We examine the efRciency of gravitational bremsstrahlung production in the process of head-on 
collision of two boosted Schwarzschild black holes. We constructed initial data for the character- 
istic initial value problem in Robinson- Trautman spacetimes, that represent two instantaneously 
stationary Schwarzschild black holes in motion towards each other with the same velocity. The 
Robinson- Trautman equation was integrated for these initial data using a numerical code based on 
the Galerkin method. The final resulting configuration is a boosted black hole with Bondi mass 
greater than the sum of the individual mass of each initial black hole. Two relevant aspects of the 
process are presented. The first relates the efficiency A of the energy extraction by gravitational 
wave emission to the mass of the final black hole. This relation is fitted by a distribution function 
of non-extensive thermostatistics with entropic parameter q ~ 1/2; the result extends and validates 
analysis based on the linearized theory of gravitational wave emission. The second is a typical 
bremsstrahlung angular pattern in the early period of emission at the wave zone, a consequence of 
the deceleration of the black holes as they coalesce; this pattern evolves to a quadrupole form for 
later times. 

PACS numbers: 



I. INTRODUCTION 



The collision of black holes, the collapse of stellar objects in the process of formation of black holes as well as the 
evolution of distorted black holes figure as promising sources of gravitational waves. The importance of these issues 
lies in the fact that the knowledge of gravitational waveforms originating from the above processes will be of crucial 
importance for the recent efforts to detect gravitational waves. However to date these issues still remain far from 
completely understood and for most situations we are forced to rely on approximation methods and numerical tech- 
niques to obtain information on wave form patterns and radiative transfer processes in the dynamics of gravitational 
wave emission 

An important astrophysical situation in which gravitational radiation is produced is that of a merger of two black 
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holes [3, in particular in a head-on collision, both with the same initial velocity. If the masses are distinct we should 
expect that the remnant of the collision would be a black hole with smaller velocity and the gravitational radiation 
produced to be typically bremsstrahlung in close analogy with electromagnetic bremsstrahlung. 

Head-on collisions of two black holes were discussed in detail, both numerically and semi-analytically, by Aninos et 
al. 0], Price and Pullin[4] and Gleiser et al.[5(]. The common feature of these papers is the use of Misner initial data@ 
to describe the initial configuration of two equal mass nonspinning black holes. Their approach bridges numerical 
relativity and perturbative techniques^ [1| to extract the gravitational wave forms at the wave zone. In Ref. [sf, I = 2 
and I = 4 waveforms at several radii of extraction were exhibited for one initial data set, corresponding to a proper 
distance between the throats for which no initial common apparent horizon is present. Also, using the extracted wave 
forms, the total gravitational wave energy output was calculated for six initial data sets corresponding to six distinct 
values of the proper distance between the throats. Initial data with small proper distance between the throats were 
also considered (initially close black holes) corresponding to a global initial apparent horizon. The latter situation of 
close black holes was examined in Ref. [J] , in which the presence of a global apparent horizon allowed for the use of 
black hole perturbation theory; the computed gravitational radiation was in accordance with the results of numerical 
computations in A comparison of the different approaches in treating the head-on collision of two black holes 
using Misner data was given in [9|] . We should also refer to Nicasio et al. |10|] as a relevant reference in the problem of 
head-on collision of two initially boosted black holes using Brill-Lindquist type datafllj. In this vein we will consider 
here the interaction of two boosted black holes moving straight toward each other along the symmetry axis with the 
same velocity, and modeled in the context of Robinson- Trautman (RT) spacetimesfl^]- In general the outcome is a 
boosted black hole with smaller velocity than the initial velocity of the two holes, and a larger amount of mass-energy 
than the sum of the individual mass-energies of the individual holes. In this process gravitational waves are emitted 
while the black holes decelerate and coalesce, producing a pattern typical of bremsstrahlung. RT dynamics is in the 
realm of characteristic initial data evolution and Bondi's mass is adopted as the total mass-energy definition. 

The plan of the paper is the following. In Section 2 the basic aspects of Robinson- Trautman spacetimes are 
presented. Section 3 is devoted to a new construction of Brill-Lindquist type initial data for Robinson- Trautman 
spacetimes which can be interpreted as two initially boosted Schwarzschild black holes; the initial data adapt properly 
to the initial value problem on null cones as is the case in RT dynamics. The main consequences of the dynamical 
evolution of these initial data are discussed in Section 4, where the numerical integration is performed using a code 
based on the Galerkin projection method. In Section 5 we summarize our main results and discuss their relevance 
and limitations as compared to previous calculations of the problem in the recent literature; an outline of future 
perspectives in this subject is also done. Throughout the paper we use units such that SttG = c — 1. 



II. ROBINSON-TRAUTMAN SPACETIMES 



Robinson- Trautman (RT) metrics are solutions of vacuum Einstein's equations representing an isolated gravitational 
radiating system. In a suitable coordinate system they have the form 

ds^ = (^X{u,e)~^-!p- + 2r^^du^ + 2dudr~r'^K'^{u,e)(de'^ +sm'^9d^^). (1) 
Einstein's vacuum equations imply 

\[u,e) = — - — + — - — cote (2) 

and 

k {\gsme)e , , 

-""^'K^^lK^^^- (3) 

In the above, a dot and a subscript 6 denote derivatives with respect to u and 6, respectively, and mg corresponds to 
the Schwarzschild mass when K = \. Throughout the paper we use units such that SttG = c = 1. The dynamics of 
the gravitational field is totally contained in the function K{u,d) and governed by Eq. ([3]), denoted RT equation. 

The initial data problem for RT spacetimes belongs to the class of characteristic initial value formulations as opposed 
to the 1-1-3 formulation, according to the classification of York[l3. The degrees of freedom of the gravitational field 
are contained in the conformal structure of parametrized 2-spheres embedded in a 3-spacelike hypersurface. For RT 
spacetimes the function K{uq, 9) given in a characteristic surface u = uq corresponds to the initial data to be evolved 
via the RT equation 
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In the semi-null local Lorentz frame given in the curvature tensor of RT spacetimes is expressed as 

„ IlABCD{u,e) IIIabcd{u,0) , Nabcd{u,9) 

HabCD = 5 1 t; \ , (4) 

where //, /// and N are objects of Petrov-type //, /// and N , respectively, displaying the peeling property [l5|. Eq. 
13] shows that RT indeed is the exterior gravitational field of a bounded configuration emitting gravitational waves, 
and for large r the spacetime looks like a gravitational wave with propagation vector d/dr. The curvature tensor 
components that contribute to the gravitational degrees of freedom transverse to the direction of propagation of the 
wave are 

i?0303 = -i?0202 = + ^ (^S) ' (5) 

where 



D{u,e) = 



Kee Ke [Kg 



(6) 



The function D contains all the information of the angular, and time dependence of the gravitational wave amplitudes 
in the wave zone. 

Now the basic equations ([2]), ([3]) have the stationary solution 

cosh 7 + cos 6 sinh 7 

where Kq and 7 are constant, and as a consequence A — Kq^ (cf. Eq. According to Bondi and Sachs[l3| this 
solution can be interpreted as a black hole boosted along the negative z-axis, with velocity parameter v = tanh7 
and mass function 771(6*) = ■moK^{9). The total mass-energy of the gravitational configuration results then in M = 
777o/2 Jq (6) sin 9 d9 = 7770 cosh 7 = ttiq/VI — v'^, where for the sake of convenience we have set Kn = 1. The 
parameter 7 is associated with the rapidity parameter of a Lorentz boost given by the K-transformations of the BMS 
groupfll]. The interpretation of ([7l) as a boosted black hole is relative to the asymptotic Lorentz frame which is the 
asymptotic rest frame of the black hole when 7 = 0. 

It will be important to exhibit a suitable expression for the total mass-energy content of the RT spacetimes or simply 
the Bondi mass. To accomplish such task it is necessary to perform a coordinate transformation to a coordinate system 
in which the metric coefficients satisfy the Bondi-Sachs boundary conditions (notice that in the RT coordinate system 
the presence of the term 2rK/K does not fulfill the appropriate boundary conditions). We basically generalize the 
procedure outlined by Foster and Newmanfrzj to treat the linearized problem, whose details can be found in Ref. 
[l8l [T9|. The result of interest is that Bondi's mass function can be written for any u as M{u,9) = moK^{u,9) + 
corrections, where the correction terms are proportional to the first and second Bondi-time derivatives of the news 
function; furthermore at the initial null surface u = uq these correction terms can be set to zero by properly eliminating 
an arbitrary function of 9 appearing in the coordinate transformations from RT coordinates to Bondi's coordinates. 
Thus the total Bondi mass of the system at u ~ uq is given by 

M(Mo) = -mo / K^{uo,e) sinOde. (8) 



III. BLACK HOLE INITIAL DATA FOR RT SPACETIMES 



The object of this Section is to construct initial data K{0, 9) for the RT equation, representing two interacting black 
holes instantaneously boosted. In order to accomplish this task we will construct 3-dim initial data for two black 
holes (that are similar to Brill-Lindquist(BL) data[ll|) from which we will extract the parametrized conformal factor 
K{9) of a family of two spheres embedded in this 3-dim geometry. We will rely on the seminal paper by MisnerQ 
which will also be a reference for our notation in the present Section. 
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Starting from bispherical coordinates [20j in the 3-dim Cartesian plane S, we are led to introduce the following 
parametrization for Cartesian coordinates 



a sin 9 sinh 77 

cosh 77 — cos 9 sinh rj 

a sin 9 sinh t] 

cosh 77 — cos 9 sinh 77 
a 

cosh 77 — cos 9 sinh 77 ' 



C0S<y9, 

sin ip^ 



z>0 



(9) 



and 



a sin sinh 77 

cosh 77 + COS 9 sinh 7; 

asin^ sinh 77 

cosh 77 + cos 9 sinh 77 
a 

cosh 77 + cos 9 sinh 77 ' 



cos if, 



sm ip, 



z <0 



(10) 



where < 77 < 00, 0<9<TT,0<ip< lix. In this parametrization, for each 77 = t/q it corresponds two spheres, one 
at ^; > and the other at z < 0. Also, for a given point P of S in z > fixed by (770, 6*0, (/^o) there is an associated 
unique antipodal point Pa — » (770, tt — Vo) in the region z < 0, so that we go from one to the other by an inversion 
through the origin (cf. Fig. 1); the Cartesian vector from the origin to the point P : [x, 7/, z) has length 



' cosh 77 + cos 9 sinh 77 



y cosh 77 — cos 9 sinh 77 
while its corresponding antipodal point is a distance r< from the origin given by 



I cosh 77 — cos 9 sinh 7/ 
cosh 77 + cos 9 sinh 77 



(11) 



(12) 



The usefulness of this parametrization will become clear in what follows. The surface z — corresponds to 7/ = 00. 
We note that the Cartesian coordinates are continuous functions, with continuous second derivatives, of (77, 0, (/?). 
Singularities occurring are the usual singularities of a spherical coordinate system. 

The flat space line element ds^ = [dx)"^ + [dyY + (dz)^ is expressed in the above parametrization as 



(cosh f) =p cos 9 sinh 77) ^ 



^77'^ + sinh^ rj{d9^ + sin^ Od^p^) 



(13) 



for z > and 2 < 0, respectively. 

We now take E as a spacelike surface of initial data, with geometry defined by the line element 



ds^ = ^'^(77 + t;o, 



drj + sinh (rj + T]o){d9 + sin 



(14) 



By assuming time-symmetric data (namely, S a maximal slice) we obtain that the Hamiltonian constraints reduce to 
— In this case the metric ([T4l) must satisfy the equation 



770)1 sinh in + Vo) 



-2||cot. + g 



sinh^(77 + 7/0 ) 

^2 



= 0, 



(15) 



where a prime denotes derivative with respect to 77. We note that the substitution K = reduces to a Laplace 
equation. Obviously the flat space solution 



FIG. 1: Projection of the spheres parametrized by rjo into the plane xz. The Cartesian vectors r> and r< localize an arbitrary 
point of the spheres = r;o in 2 > and z < 0, respectively. 



K{7^ + 770, = ^ , , . , . ^ , (16) 

cosh(77 + 770) T cos y smh(?7 + rjo) 

satisfies Eq. (fT5| . from which it follows that the function 

K{V, 9)=( , + \ (17) 

y ^cosh(r7 + 7;o) — cos6'sinh(j7 + 7?o) \/ cosh(?7 + r;o) + cos sinh(7? + r;o) j 

is a nonflat solution of (|15p . where ai, a2 and ryo are arbitrary positive constants. We may interpret the non-flat 
3-dim metric defined by ITTI) . 

ds^ = a (j],e) d-q'^ + smh^ {7^ + 'no){de^ + sin edip'^) , (18) 

as a BL-type solution given in bispherical coordinates. In fact a straightforward manipulation shows that the metric 
PS)) can be rewritten as 

4 



2 1/ , /cosh(r7 + 770) + cos6'sinh(77 + ryo) \ ,2 

ds = - a2+aiW TT — ^ { . . ' — ■ 't dSflat{+) 

2 \ y cosh(r; + 770) — cos t/ sinh(77 + ryo) / 



1 / /cosh(77 + r;o) - cos6'sinh(j7 + 7;o) \ , 2 ,,rix 

2 \^ y cosh(j7 + r;o) + cos 6'smh(j7 + T^o) / 

where c?Sflat(+) ^'^'^ '^•®flat(-) respectively the form of the flat space metric for the z > and z < domains. For 
r] » rjQ and corresponding x,y, z » a we may express 

ds^ = i a2 + ^ ) d4at(+) + i ( ai + ^ ) dsflat(-) , (20) 
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and returning to Cartesian coordinates the 3-geometry can be given in the approximate form 

This allows us to interpret the initial data (|17p as two interacting Schwarzschild black holes instantaneously at 
rest, localized at ?'>(r/o) and ?'<(?7o), in z > and z < 0, with masses Mi =^ 2|q;i'^ 0:2/(0:1^ + a2'*)| and M2 — 

2|a2^ai/(ai^ + ck2^)| in units of a, respectively. It is worth noticing that the radial Schwarzschild isotropic type 

coordinates r>(?7o) and ry{r]o) are functions of the bispherical coordinate 9, and determine the Euclidean distance 
of points of the spheres rj — rjQ to the origin. The minimal Euclidean distance between the two spheres is given by 
L = (r> + r<) evaluated along the z-axis, resulting in L = 2a exp(— r^o)- 

From the above construction we can now extract initial data for the RT dynamics, which has its initial value problem 
on null cones. We note that the geometry of the two spheres located at rj = i]q (cf. Fig. 1) contains all the information 
on the initial data for vacuum field equations through the conformal function (jl7p calculated at 77 = 0. Based on the 
initial data formulation on characteristic surfaces proposed by D'Inverno and Stacheljlil [22|i in which the degrees 
of freedom of the vacuum gravitational field are contained in the conformal structure of 2-spheres embedded in a 
3-spacelike surface, we are then led to adopt the conformal structure given by 

K{7jo,e)^( , + ^ (22) 

V V cosh r;o — cos Q sinh 770 v cosh 770 + cos Q sinh 770 / 

as initial data for two interacting Schwarzschild black holes to be extended along null bicharacteristics and propagated 
along a timelike congruence of the spacetime. A restricted spacetime of two interacting Schwarzschild black holes 
may then be constructed locally as the product of the two 2-sphere geometry times a timelike plane (u, f) generated 
by a null vector djdf and a timelike vector djdu with geometry da^ = c?(u^ f, ff) dv? + 2 dr. The four geometry 
of the product space is then taken as 

da' = Q^(ii, f , 6^) d%i' + 2 du dr — K^^rjo, 6, u)(d6'^ + sin^ Odifi^). (23) 



Eq. I|23p is the Robinson- Trautman metric, the dynamics of which (ruled by Einstein's vacuum field equations) 
corresponds to propagating the initial data K{r]o, 9,u = 0) = K{riQ, 0) (cf. (I22p ) forward in time from the characteristic 
initial surface u = uq. 

This characteristic initial data problem presents an ambiguity in the interpretation of the parameter r/o- On one 
hand it may be interpreted as the instantaneous boost parameter of the black holes along the z axis; on the other 
hand it may work as defining a distance between the two black holes. A reason for this is that the family of 2-spheres 
parametrized with 770 (via the conformal factor (|22|)) are propagated along the timclike vector d/du that is boosted 
relative to the timelike vector d/dt. For the simple cases of ai = or q;2 = this boost transformation is generated 
by the conformal factor (|22p itself [1^. We therefore are led to interpret (|^^ as representing two initially boosted 
black holes, with opposite velocities v = tanhr7o along the z axis and with an initial Euclidean distance parameter 
L/a = 2 cxp(— 770). 



IV. NUMERICAL RESULTS: EMISSION OF GRAVITATIONAL WAVES AND MASS LOSS 

We now evolve the initial data (1221) via the RT equation from u = Q. This equation is integrated numerically using 
the Galerkin method (see Refs. [14|.[2^ for details) in which the approximate solution of the RT equation has the 
following form 

K{u,e) = Ao exp l^i^ 6fe(«)P(cose)^ (24) 

where Aq is an arbitrary constant, hk{u) are unknown modal coefficients and Pk{cos0) are Legendre polynomials. 
According to the Galerkin method a set of -|- 1 ordinary differential equations determines the evolution of the modal 
coefficients, and therefore the function K{u,6). The initial conditions bk{Q), k — 0,1,. .,N must correspond to the 
initial data (|^^ and are evaluated from 
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_ (2 In {K{x)Ao^),Pk{cose)) 
{Pk{cos9),Pk{cos9)) 



(25) 



where the brackets define the orthogonal product in the projection space of Legendre polynomial normalized according 
to {Pk{cos 9), Pj {cos 9)) = Pk{cos9),Pj{cos9)sm9d9 = 2Skj/i2k+l). The effect of increasing the truncation order 
TV is illustrated in Fig. 2 with the distribution of the absolute error between the exact and approximate initial data, 
respectively given by ((221) and i^approx = exp J2k=o ^fc(0)P(cos 6*)^ . 



.1e-1 
.1e-2 : 
.1e-3 : 

g 

t*(i 1e-06 

i 

^ 1e-07 

1e-08 1 
1e-09i 




-0.5 



0.5 



FIG. 2: Plot of the absolute error between the approximate and exact initial data for a2 = 1.0, ai = 0.2 and 7 = 0.6. The 
curves correspond to truncations orders A'' = 4, 7, 11. Here we have replaced the angular coordinated 6 hy x — cos 9. 

We have done exhaustive numerical experiments considering truncation order iV = 13 and taking into account 
distinct values for ai, a2 and 779. Basically, the numerical outcomes can be summarized as: 

• ai ^ a2 (two unequal mass black holes): the final configuration is a boosted black hole with a smaller velocity, 
or r]o final < ?7o; 

• ai — 0,2 (two equal mass black holes): the final configuration is a black hole at rest, or rjo final = 0. 
The final configuration iffinai can be reconstruct as 



if final (S) = exp i ^ bk{ut)Pu{c 



(26) 



where at the final time of integration u — Uf, the modal coefficients bk{u{) ~ constant (at least up to 10 ). 

Kq final 



Kfin.ii9) 



cosh r]o final + COS 9 Sinh l]o fina 



(27) 



It is worth noting that for all if gnai it is always possible to find values for Kq gnai and 770 final such as the inequality 

\Kfina.l{9) — ifo final/ (cosh 7?o final + COS 6* sinh 7?0final)| < 10"^ (28) 

holds for all < < TT. Therefore the final configuration of the system is a boosted black hole; the final velocity and 
rest mass are determined and given respectively by Ufinai = tanh(77o final) and rrifinai = "^oifofinai- 

One of the most interesting aspects of the dynamics of the RT spacetimes is the emission of gravitational waves 
and consequently the mass loss of the initial configuration. Let us explore this feature considering the situation in 
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which a2 and 770 are fixed and ai is the free parameter. Therefore at u = there are two boosted black holes directed 
towards each other with the same velocity and separated by a fixed distance. One of them has fixed mass while the 
mass of the second varies by changing ai. In the numerical experiments we have considered the range for which 
1 <^ a^joLx « 1. In particular for a^joix <^ 1 the initial data can be recast as 



cosh rio + cos 6 sinh rjo (cosh^ r]o — cos^ 6 sinh^ ??o)^/^ ' 

meaning that it can be viewed as a boosted black hole with total mass-energy E2 = moaf cosh 779 perturbed by a 
smaller black hole whose total mass-energy is Ei — moaf cosh 770 directed towards the first black hole. Noticed that 
since Ei <C E2 its contribution to the total mass-energy associated to the initial data can be neglected. 

According to Eq. ([5]) the total mass-energy of the system evaluated at u — is given by Eq = 
1 /2r7io Jq (^) sin 9d6, so that the efficiency of the process of emission of gravitational radiation can defined bv[27| 

^^ Eo-MsH ^ (30) 



where Mbh is the total mass-energy of the resulting final black hole. In the present context part of the kinetic energy 
and the interacting energy of both black holes is radiated away, and another part is absorbed by the first black hole 
increasing its rest mass to the amount "loi^o final- ^he RT equation is integrated taking into account the initial data 
([22|) with 770 = 0.3, a2 — 1 and varying ai inside the range [10"'^, 1.3]. For each value of ai the final outcome is 
determined and identified as a boosted black hole with smaller velocity, Ufinai = tanh 770 fjjjai < tanh 770 and total 
mass-energy such that Mbh > Ei + E2 as expected from the area theorem of black holes[25|. In Fig. 3 the result 
of the numerical experiments for which 10"'^ < ai < 0.3 is displayed by the plot of the efficiency A as function of 
Mbh, along with a continuous line accounting for an analytical function that fits quite well the numerically generated 
points. Then, as we have shown previously [la] the function A = A{Mbh) inspired from non-extensive statistics [26l] 
is given by 



A = A,nax (1 - y^)'/'^^ (31) 

where y = E2/MBH is the ratio between the initial mass-energy of the first black hole (fixed parameter a2) and 
the mass-energy of the resulting black hole, A^ax is the maximum efficiency of the process attained in the limit 
Mbh 3> 7 and q are the free parameters that characterize the non-extensive ;26i] relation. The best fit shown in 
Fig. 3 demands 7 ~ 0.525, q ~ 0.502 and A^ax — 0.000977, which is a quite small efficiency even for a considerably 
high initial boost of tanh(0.3) « 0.3 or about 90,000 km/s. As a matter of fact we expect that the maximum efficiency 
will depend on the initial boost of both black holes. In any case the maximum efficiency is in accordance with the 
value about 0.07% found by Smarr and Eppleyfl^ for the collision of two equal and non-rotating black holes in a more 
general spacetime. Although the non-extensive relation (|3ip is in excellent agreement with the numerical results in the 
range of ai under consideration that covers the interval < Mbh ^ 5.0 M^^, where Mg^ is the minimum value 

of the mass of the final black hole, it fails to describe the entire set of points. The numerical value of the maximum 
efficiency, 0.000598, obtained when ai = a2 is distinct from the one predicted by Eq. ([3T|) . Also for ai > a2 the 
efficiency seems to decrease, as would be expected from the symmetry in the interchange of ai and a2 in the initial 
dafe ISMSidlfM limitation of the non-extensive relation over entire range of numerical points, we cannot underestimate 
the fact expressed by the fitting of Fig. 3. Let us briefly discuss some implications of relation ([3T|) . starting from the 
domain of A/Amax ^ 1, or relatively very small efficiency. In this situation it is clear that E2 ^ Mbh and therefore 
y ^ 1, so that Eq. (pij) can be written as the following scaling relation 



Based on the linearized theory of gravitational wave emission we can accordingly demonstrate that the parameter q 
must be 1/2, the same value obtained from our numerical experiments in the full nonlinear regime of RT dynamics. 
Actually, in the limit A/Amax <C 1 we might consider the initial mass given by Eq — E2 + SE, where 6E/E2 = e <C 1. 
From the quadrupole formulafSsJ the mass loss can be written as d6E/dt ^ 0{e^) which is negative definite. This 
means that the perturbation falls into the hole such that at the end its mass becomes Mbh — E2 + (^-Eabs, where 
(5i?abs is the total amount of mass absorbed. The other fraction of the perturbation (5i?rad is radiated away, and as a 
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FIG. 3: Plot of the ratio Mbh /ma versus the efficiency A in the process of gravitational wave emission due to the collision of 
two black holes for 10~^ < cti < 0.3 and rjo — 0.3. The continuous line is the best fit of the points by the function (|31[) with 
q ~ 1/2. 



consequence of the quadrupole formula it follows that <5i?rad/-E'2 ^ 0{e^). Taking into account these relations we can 
express the efficiency as A = {SE — 5Eahs)lE2 ^ 0{e^), and — 1 ^ C>{e)^ which leads to 



i - 1 oc Ai/2 
V 



Mbh = -52(1 + const. 



(33) 



Comparing Eqs. ([32|) and ([33|) we find q = 1/2. It must be emphasized that the above derivation is general and 
therefore not restricted to the realm of Robinson- Trautman spacetimes. 





-0.6 -0.4 -D.2 



0.2 0.4 0.6 0.8 



-0.008 -0.006 -0.004 -0.002 0.002 0.004 0.006 0.008 



FIG. 4: Initial amplitude of gravitational waves \D[u = 0, a;)| for which a\ — 0.3, 02 — 1 and rjo = 0.6. The second figure shows 
a sequence of polar plots of \D{u,x)\ at early times. Notice that the angular pattern is typical of bremsstrahlung analogous 
to the electromagnetic bremsstrahlung of a decelerated charge along its direction of motion, where the lobes open due to the 
deceleration. For later times the angular pattern evolves to a typical quadrupole form. 



The second aspect we will discuss is the wave zone angular pattern of the gravitational waves emitted as described 
by the function D{u,9) (cf. Eq.Q). Its form is valid for any value of r sufficiently large since 1/r is a multiplicative 
factor to D{u, 9) in the wave zone curvature component. As we have mentioned the case ai = a2 does not produce a 
boosted black hole due to the symmetry of the initial data. In order to illustrate the type of angular pattern of emitted 
waves we set 770 — 0.6 corresponding a considerably high velocity and two black holes with initial masses of about the 
same order, or 02 = 1, ai = 0.3. In Fig. 4(a) we show the initial pattern of these two black holes represented by the 
plot of \D(u = 0, a;)| versus x = cos 0, that exhibits a dominant emission in the northern hemisphere, with a maximum 



10 



amplitude at a; ~ 0.81. We note that the emission cone of maximum ampUtude opens up due to the decrease in 
the deceleration of the system; also as it opens up the maximum amplitude decreases and the pattern evolves to a 
quadrupole form for later times. The final state is that of a boosted Schwarzschild black hole with rio fi„ai ~ 0.355. 

V. CONCLUSIONS AND FINAL DISCUSSIONS 

In this paper we have studied the head-on collision of two non-spinning black holes in the realm of Robinson- 
Trautman spacetimes. Brill-Lindquist type initial data were constructed that represent instantaneously two 
Schwarzschild black holes directed towards each other with the same velocity. The RT equation was integrated 
numerically using a combination of Galerkin and collocation methods [24] for black holes with distinct initial masses. 
In general the resulting configuration is a boosted black hole as described by Eq. ([27|) characterized by a smaller veloc- 
ity than both initial velocities, while the final mass is greater than the sum of the initial masses, or Mbh > Ei + E2- In 
this process a fraction of the total initial mass is extracted by gravitational waves, whose amount has been evaluated 
for distinct values of the parameters. 

Besides the derivation of the initial data the paper presents two relevant aspects of the dynamics in the full nonlinear 
regime. The first is a relation between the efficiency A (cf. Eq. [3T|) of the gravitational wave extraction and the mass 
of the final black hole. The remarkable feature is that the numerical results expressed by the points {Mbh, A) could 
be fitted by an analytical function inspired from the non-extensive statistics as we have done considering other initial 
data, with entropic index q ~ 1/2. The question of whether or not such a relation is an artifact of Robinson- Trautman 
spacetimes is crucial for guiding us to further investigation on this issue. Nonetheless we have guaranteed that at 
least for very small efficiency the power-law p2p derived from the quadrupole formula is general (not restricted to 
Robinson- Trautman spacetimes) and appears as a well-defined limit of the non-extensive relation (|3ip with q — 1/2. 
Another important aspect we have displayed is the bremsstrahlung pattern of gravitational radiation in the wave 
zone, a consequence of the deceleration of both black holes as they coalesce. 

As compared to previous works (Refs. 0] [1| Q [1] most of which deals with the head-on collision of two 
Schwarzschild black holes, our approach differs basically in that we have adopted the characteristic surface initial data 
formalism. The latter indeed has several advantages for the description of the gravitational radiation and also for 
the construction of marching algorithms [^oj . Therefore, we consider that the problem of collision of two black holes 
in the realm of characteristics is worth studying. In this direction, an accurate code based on the Galerkin method 
with collocation was constructed to integrate the field equations. The code is highly stable for long time runs in the 
full nonlinear regime so that we are able to reach numerically (up to 10"'', cf. Eq. (j28p ) the final configuration of 
the system, when the gravitational emission ceases. The final configuration is a remnant Schwarzschild boosted black 
hole, with smaller velocity parameter (jyofinai < Vo), and definite rest mass mfinai — "^o^ofinai ^^'^ total mass-energy 
Mbh — (^™o-?^o final cosh77o final j (cf. Eq. (US])). We were therefore able to evaluate the efficiency of gravitational 
radiation emission in extracting mass-energy of the source in the full nonlinear regime as a function of the final black 
hole mass and the initial mass-energy of the system. This procedure is in the line of, and consistent with the general 
estimates made by Eardley[27]. Our approach is alternative to the one made in the above mentioned references 
on head-on black holes collisions, which evaluated the asymptotic energy flux carried off by the gravitational wave 
(or equivalently the efficiency of the gravitational wave extraction of mass from the source) through the use of the 
extracted gauge invariant Zerilli function[9j] . The comparison of the results in the two procedures is not straightforward 
due to the fact that the one parameter Misner data and our three-parameter BL-type data are somewhat distinct, and 
numerical evaluations of the efficiency were made with variation of distinct parameters (for instance initial separation 
or mass-energy of the final configuration). We should mention that a few numerical tests made with fixed ai — 1.0 
and a2 — 0.1 and varying the boost parameter 770 in our approach showed a saturation tendency of the efficiency for 
increasing ryo, which is in accordance qualitatively with the second order perturbation results of Ref. [l0[. We expect 
to examine these issues in the future. 

Finally, we may cite further steps in our future investigations: (i) extension of the current work by constructing and 
evolving initial data that represent a non-central collision of two black holes; (ii) the construction of characteristic 
initial data that may accommodate spinning black holes; (iii) a more general and detailed study of the amount of mass 
extracted by gravitational waves and (iv) the generality of the non-extensive relation pTjl connecting the efficiency 
and the mass of the final black hole. 
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